•CLASSIFIED 


T"~ ~W~ 



I 


f 


Rsf2'/yi40§d 

Iu< itie 
§ 

{VICES TECliNICU. INFORMATION AffiNCY 
1UNGM H.iLL STATION 
Z IWJK i;:, VIRGINIA 




^CLASSIFIED 









NOTICE: When government or otheraiitravings, speci¬ 
fications or other data are used, ifer any purpose 
other than in connection with a 4M±nltely related, 
government procurement operation, ^ tQae U. S. 

Government thereby incurs no respcijtcm.Bibility, nor any 
obligation whatsoever; and. the fafKctt. that the Govern¬ 
ment nay have fonulated, furnishes!- sc,, or in any way 
supplied the said drawings, spec±J»ffications, or other 
data is not to he regarded by impa;ilcation or other¬ 
wise as in any manner Ucensing'haHne holder or any 
other person or corporation,/or aammreying any rights 
or permission to manufacture, userai ior sell any 
patented invention that'may in araifij- way be related 
thereto. ^/ 



/ 


o 

\Q 


MASSACHUSETTS INSTITUTE OF TECHNOLOGY 
NAVAL SUPERSONIC LABORATORY 



Technical Report 387 

INFRARED RADIATION FROM 
AERODYNAMICALLY HEATED SURFACES 


Lawrence D. Lorah 


June 1959 


A S T f A 

fiBgBafPfBflraa 

HIV «• tsMi ' 

F« 

i. uisrarerrsaJ 


"*• c* 






op 


(s !" 3 l / 

XEROX 




Massachusetts Institute of Technology 
Naval Supersonic Laboratory 
Cambridge, Massachusetts 


Technical Report 387 


INFRARED RADIATION FROM 
AERODYNAMICALLY HEATED SURFACES 

by 

Lawrence D. Lorah 

Contract No. Nonr 1841(40) 

MIT DSR No. 7683 


June 1959 




ABSTRACT 


A theoretical method is presented for determining a 
spectral and spatial description of the infrared radiation from 
aerodynamically heated surfaces. The basic laws of infrared 
emission from opaque surfaces are reviewed and put into a 
form useful in this application. Heat transfer rates to flat 
plates and cones in laminar and turbulent supersonic flow are 
examined and pertinent equations are pre sented. Solutions for 
equilibrium temperatures are obtained for representative 
flight conditions and the radiation laws are applied to produce 
characteristic radiation patterns at various altitudes and Mach 
numbers. As an example, the radiation from the skin of a 
hypothetical aircraft is found for various flight conditions. 
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SECTION I 

INTRODUCTION 

As infrared sensitive devices become of greater and greater 
importance in aerial warfare, it is imperative that the equipment 
designer have,at hand,knowledge of target radiation signatures. Even 
though such considerations as background characteristics and optical 
design help determine the ultimate capabilities of a detecting or tracking 
system, one of the most basic pieces of information necessary for 
intelligent design is the level and spectral distribution of the infrared 
energy coming from the target. In the case of aerial craft this infrared 
radiation can come from several sources. The structure becomes 
heated and the exposed surfaces emit radiation peculiar to their temper¬ 
ature distributions and surface characteristics. The propulsion system 
produces a strongly radiating plume of gas if it is a jet or rocket 
system. Sunlight and earth radiation can be reflected from the surfaces, 
and at very high speeds and at extreme altitudes the heated air layer 
surrounding the craft and the hot wake may produce significant infrared 
energy. 

The purpose of this paper is to examine the first source of 
radiation, the aerodynamically heated, external surfaces of the vehicle, 
and develop satisfactory methods of estimating its level and spectral 
distribution. It has proven difficult to obtain satisfactory radiation 
measurements of heated airframes, especially under realistic flight 
conditions. Even if methods were developed to make accurate measure¬ 
ments on the pertinent vehicles, a potential enemy would not allow measure 
ments to be made on his weapons. Therefore it is necessary to develop 
engineering methods of predicting the radiation characteristics of these 
surfaces once adequate target models and operating conditions are 
established. High accuracy in the calculations is not usually necessary 
due to the uncertainty in target geometry and operating conditions. 
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However, care should be taken so that the results obtained are meaning¬ 
ful and should be as accurate as the target characteristic hypothesis. 

The problem is approached by first reviewing the basic laws of 
infrared emission from opaque surfaces in order to clarify nomenclature 
and to establish useful equations and approximations. Then the heat 
transfer rates and equilibrium temperatures for the surfaces under con¬ 
sideration are examined. The major mode of the heat transfer at high 
speeds at reasonable altitudes is aerodynamic heating and it is investigated 
here in the continuum flow region. The heating analysis is restricted to 
simple shapes (cones and flat plates) from which model targets can be 
constructed. The effects of radiation heat loss as well as heat input by 
solar and earth radiation are discussed. No account is made of possible 
internal heat sources or sinks. The laws of radiation are then coupled 
with the temperature analysis to produce radiation patterns at various 
conditions. Several approximations are then evaluated and general trends 
in radiation characteristics with varying conditions are discussed. The 
entire discussion is restricted to radiation at the target,and the effects 
of atmospheric attenuation are not considered. 
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SECTION II 

EMISSION OF INFRARED RADIATION 

In beginning this discussion of infrared radiation from heated 
surfaces it is important that the basic laws and concepts of thermal 
emission be reviewed. The purposes of this section are, to establish 
a radiation nomenclature to be used throughout the paper, to establish 
the basic equations in a form convenient for later use, and to make 
clear the necessary assumptions and approximations used. The 
following is a very brief discussion of the spectral and spatial distri¬ 
bution of infrared radiation from an opaque surface as a function of 
surface characteristics and temperature. For a more complete and 
comprehensive treatment the reader is referred to standard physics 
or heat transfer texts such as Refs. 1 and 2. 

In discussing infrared radiation use will be made of a " black 
body" which is perfectly absorbing at all wavelengths-one which reflects 
or transmits no radiation. As an example, a small aperture in a 
large cavity in thermal equilibrium satisfies the requirement of perfect 
absorption, and is thus a "black" surface. For such a surface, Stefan's 
law states that the total flux (power per unit area) radiated is proportional 
to the fourth power of the absolute temperature of the surface that is. 


l bb 


= <r T‘ 


( 2 . 1 ) 


where 

q^b = total flux from black body 

T = temperature 

c = Stefan-Boltzman constant 
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A black body is not a good representation of actual surfaces, and 
to account for this departure from ideal conditions the ratio of the total 
flux radiated from a real surface at a given temperature, q(T), to 
the total flux radiated from a black body at the same temperature is 
defined as the emissivity, « (T); 


-9-UJ- = e (T) 

%b (T > 

or using (2.1) 

q(T) = c(T) a T 4 


( 2 . 2 ) 


This temperature dependent, total hemispherical emissivity is defined 
using q(T) and q bb ( T ) evaluated over all wavelengths, and it is assumed 
that c is constant over the entire spectrum. A more rigorous treatment 
would include the variation of e with wavelength, a variation which becomes 
significant when considering gaseous radiation, or radiation from partially 
transmitting solids. However, this is of less importance in dealing with 
opaque surfaces of structural material, since the variation of c with 
wavelength is small and can be neglected in most cases. All surfaces in the 
following discussions will be considered "grey" (i.e. , c not a function of 
wavelength), but the strong variation of emissivity with temperature must 
be accounted for. This is probably best done by referring to measurements 
of emissivity on structural materials made in the temperature region of 
interest. 

Lambert's cosine law describes the spatial distribution of the 
radiation from an ideally grey surface element. Referring to the 

NORMAL 
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diagram, Lambert's law states that the increment of power, dw, 
radiating from the element of area dA at temperature, T, into the 
element of solid angle, df2, can be written, 


dw = 


« <r T 4 


cos 8 d£2 dA, 


(2-3) 


where 9 is the angle between the surface normal and the direction of 
dw. Upon integrating (2.3) over an entire hemisphere, (2.2) is obtained 
since 


dw 

dA = q • 

The situation of interest here is one of collecting the surface radiation 
at long ranges with relatively small apertures. Therefore the solid angle 
under consideration is very small, so cos0 does not vary significantly 
over and will be assumed constant for any given surface element aspect. 
Using this assumption, (2. 3) is rewritten in the more convenient form 
of flux (power per unit of radiating surface area) per steradian subtended 
by the collecting aperture; 


dq e (r T 4 

dfi = -7- COS0 (2.4) 

Now, (2.4) holds only for a point source, but the assumption of target 
dimensions much less than the detection range allows the complete 
target to be considered a point source. The total power of the target 
as a point source is found by multiplying (2.4) by dA and integrating over 
the target surface area in view. 

“ST = J Tfl- dA = ~~ j T4 cose dA (2.5) 

It should be emphasized at this point that Lambert's law as shown 
in (2.3) and (2.4) holds only fd'r ideally grey surfaces on which e is 
constant for all observation angles. For real surfaces there is a 
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disparity in the spatial distribution of radiant energy as the surface is 
observed from angles near 90° to the surface normal. This variation 
is best described by allowing the local emissivity to be a function of 6. 

In dielectric materials this emissivity falls off to a theoretical value of 
0 at 90° to the normal. Conductors display an opposite variation, 
theoretically the emissivity tends to increase as 6 approaches 90 . The 
diagram indicates the nature of the variations by a polar plot of e (0) 

€( 0 ) 



versus 0. These departures from ideal conditions are discussed in some 
detail in Ref. 2 . However, they will be neglected in this paper for the 
following reasons. First, the effects are small when the emissivity is 
relatively high, as it is in normal structural materials • And, second, 
the correction calculations cannot be made with any degree of accuracy 
anyway because of the uncertainty in the surface characteristics of a 
target (smoothness, oxidation, etc. ) and the lack of data on the physical 
characteristics of such surfaces even if they were described (index of 
refraction for long wavelengths, conductivity, permeability, etc.). Some 
discretion must be exercised when the line of sight is nearly perpendic¬ 
ular to the surface normal (greater then 75° or 80°), but it can be seen 
that the results of neglecting the real surface effects tend to be conserv¬ 
ative, at least in the case of metallic surfaces. 

Up to this point the flux integrated over wavelengths from 0 to oo 
has been considered, however, it is important to know the spectral 
distribution of the target radiation in order to select sensitive elements 
and account for atmospheric absorption. The most convenient, continuous 
description of energy as a function of wavelength utilizes the derivative 
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of energy with respect to wavelength. If the flux emitted in the wave 
length interval between X and X + dX is defined as q x d\, (or more 
simply, q x = dq / d\), Planck- s law gives the relation; 


2uhc* * -- dX . (2.6) 

q Xbb ( x .T) dx = x5 ^ch/XkT^ 

This is the basic equation giving q x for a black body as a function of 
wavelength and temperature. Integrating (2.6) from X - 0 to X - oo 
yields (2. 1) with, 

2 Tt^k 4 

= 15 c* h~ * 

In order to obtain a complete description of the radiation from 
surfaces as a function of temperature, wavelength and geometry (assuming 
long ranges, relatively small apertures, and ideally grey surfaces) 
Lambert- s cosine law and (2. 2) are incorporated into (2.6) to yield, 


dq X 

d n 


X 5 


2 h c z t _ 

/ch/kXT_ ^ 


cos 0 


(2.7) 


Multiplication of (2.7) by dX and then integration between two wavelengths 
gives the flux emitted between those two wavelengths per unit solid angle. 
The form of (2. 7) is convenient in that the energy of various targets can 
be compared directly by evaluation of the right hand side which is independ¬ 
ent of range and aperture size- The latter parameters are easily dealt 

■with since, 


n = A / R 2 . 

c' 

The appendix gives a simple way of producing plots of (2.6) or (2.7) for 
any temperature. 
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If an analytic expression for - is desired, a great reduction 

in complexity can be obtained be approximating Planck's equation, (2.7), 
by Wien's formula. This can be derived from (2.7) by assuming that 
ch/ X.Tk » 1 so that the -1 in the denominator can be neglected. This 
yields 

* « 2h c 2 «- ch / XkT cos 6 (2-8) 

dw 

This approximation holds well when X.T « 1.4 cm'K. 

Now consider a surface radiantly exchanging heat with its 
surroundings. According to Kirchhoff' s law the fraction of incident 
energy absorbed by an opaque surface is equal to its emissivity. So, 
denoting the surface under study and the external radiation sources by 
the subscripts 1 and 2 respectively, the net radiation heat flux from the 
surface can be written, 

[ T| . -W,v] ,2 ' 9) 

The k's are geometrical factors accounting for the extent and orientation 
of the external radiation sources. The strongest external sources of 
radiation in the present heating problem are the atmosphere and the sun. 
There is, of course, a certain amount of " earthshine"; however, it will be 
assumed that the intervening atmosphere has absorbed enough of this 
radiation so that it can be neglected. It has been shown (Ref. 3) that 
optically thick layers of clear air as well as clouds and haze radiate 
essentially as a black surface at the ambient air temperature. For air 
totally surrounding the surface under consideration, the geometrical 
factor equals one. The solar energy can be accounted for by using the 
common solar constant and a trigonometric factor to account for the 
surface's orientation. Explicitly, 

q 8 = ocosB 

a = 0. 1395 watts/ cm 2 


TR 387 


8 


Now (2.9) can be written to include these two sources of heat; 


q = <r «. ["t 4 - T 4 - cos9 1 (2. 10) 

M ill air i t sj 

Obviously the last two terms in (2. 10) can be neglected if the temperature 
of the surface is high enough. 
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SECTION III 


CONTINUUM FLOW AERODYNAMIC HEAT TRANSFER 

In this section the equilibrium temperatures and heat transfer 
rates of aerodynamically heated surfaces in continuum flow will be in¬ 
vestigated. This flow regime embraces the region of " convention.'; 1" 
flight where the air can be considered an isotropic continuum and the 
long developed laws of fluid dynamics can be applied. The temperatures 
and heating rates on two simple shapes from which most model targets 
can be constructed (flat plates and cones) are first considered by rela¬ 
tively sophisticated methods to obtain temperature distributions as a 
function of distance from the leading edge. Then, less laborious, ap¬ 
proximate solutions are presented to obtain an isothermal approximation 
to the temperature distribution. Since aerodynamic heat transfer is a 
boundary layer phenomena, the two classes of boundary layer flow are 
considered, laminar and turbulent. 

In the following discussion, when local free stream conditions 
are referred to, it is to be understood that these are the conditions just 
outside the boundary layer at the point in question. For the simplest 
case - the flat plate at zero angle of attack - these are the undisturbed 
ambient conditions. However, in the case of a flat plate at angles of 
attack other than zero and in the case of a cone, account must be taken 
of the changes in temperature, pressure and velocity as the flow passes 
through a shock wave or an expansion wave. For attached -nock con¬ 
ditions, wedge and conical shock solutions such as those appearing in 
Ref. 4 can be used to find the local Mach number and pressure behind 
these discontinuities. The fact that the stagnation temperature remains 
constant through an uncurved shock wave yields local static temperature 
through the Mach num ber - temperature relation (again found in Ref. 4). 
For surfaces at a negative angle of attack, the Prandtl-Meyer expansion 
method serves the same purpose, and is briefly discussed in the same 
reference. 

To carry out the calculations in this section some of the physical 
properties of air must be known. The ambient conditions as a function 


of altitude must be known, and the properties (viscosity, conductivity, 
specific heat at constant pressure, etc. ) as a function of temperature 
and pressure are also required A useful source of this information 
is Ref. 5. 

A quantity which tends to appear in most continuum flow heat 
transfer calculations is the adiabatic wall temperature, T aw . This is 
the wall temperature at which there is no aerodynamic heat transfer 
between the wall and the stream. For laminar flow it has been found 
that, 

T aw i =T ° [ l + ~ X T~ Prl/ * M *] ’ (3-D 

and in turbulent flow, the expression becomes 

T aw t = T “ [ 1+Jc r- Prl/s MZ ] • (3.2) 

The adiabatic wall temperature is often found through the use of a 
recovery factor, r, which is defined 

(3.3) 

The recovery temperature, or adiabatic wall temperature, is the same 
in either system with r = Pr^ 2 in laminar flow and r = Pr 1 ^ 3 in turbulent 
flow. 

Laminar Flow Over a Flat Plate 

The first analysis of the laminar flat plate used in this section 
is that developed by Lighthill in Ref. 6. In this analysis the solution for 
compressible, laminar boundary layer with heat transfer and zero pressure 
gradient found by Pohlhausen is used. The original solution is restricted 
to the case of a uniform wall temperature; however, Lighthill modified 
it to apply to a non-uniform wall temperature by adding a solution to the 
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heat conduction equation in the boundary layer. The resulting expression 
for the aerodynamic heat transfer as a function of the distance from the 
leading edge, x, is, 

q (x) = 0.339 K Pr‘/ J ( (- x ^ ) 


r 


(X 






(3.4) 


where the subscript 1 denotes the variable of integration. Expression 
(3.4) is set equal to the total radiation heat loss, neglecting back 
radiation from the surroundings (2. 2) and after a great deal of non¬ 
dimens ionalization and application of boundary conditions, the solution 
to (3.4) is shown to take the form, 

T(x) =T aw [F(Z)] . (3-5 

In this expression, 

- - (»* (+-)*• 

and the reference temperature, T v , is given by 


T* ’ . 0-678 K_ / ui .j , (3 .7) 

€ (T K \ P- ) 

The function F ( Z) in (3.5) has been evaluated, and is given by the 
infinite series, 

F(Z) =1 - 1.461 Z + 7.252 7 - 46.46 Z + 332.9 Z 4 -2538 Z* 

6 

+ 20120 Z -- (3.8a) 

for Z < • 1, and for Z > . 5, 


F(Z) =0.8409 Z 


-V4_ 


0.1524 Z 


-Vz _ 


0. 0195 Z 


- 3 / 4 _ 


0.0038 Z 


— I 


( 3.8b) 


The portion . 1 < Z< . 5 is simply a curve faired between the two 
solutions since no singularity should exist in that region. The function 
F ( Z) versus Z is shown in Fig. 1. In evaluating the above equations 
the adiabatic wall temperature, T aw , and the local velocity, u, are 
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computed for local free stream conditions just outside the boundary 
layer. The air Prandtl number, Pr, thermal conductivity, K, 
viscosity, p,, and density, p, are evaluated at the wall temperature. 

The pressure is considered constant throughout the boundary layer. The 
analysis can be extended to flat plates at angle of attack since there is 
no pressure gradient over such plates in supersonic flow. 

By returning to Pohlhausen's original isothermal wall solution, 
an approximate temperature for the plate can be obtained. The engi¬ 
neering approach of a heat transfer coefficient, h, and a temperature 
difference, or temperature potential, is used to estimate the aero- 
dynamically produced heat flux, i. e. 


q = h (T - T) 
n aw 


(3.9) 


For the laminar case the adiabatic wall temperature is given by (3. 1), 
and Pohlhausen, as quoted in Ref. 6, gives the local heat transfer 
coefficient as a function of the streamwise distance from the leading 
edge, x, as, 


h = 0. 332 Pr /3 K ( 


7s 




(J.X 


)'/z 


(3. 10) 


Again the physical properties of conductivity, viscosity and density are 
evaluated at the wall temperature; the velocity and adiabatic wall temper¬ 
ature are evaluated at local free stream conditions. 

To obtain an average heat input over the length of the plate,(3.9) 
is integrated over the plate length and divided by the total plate area. 
Specifically, considering a plate of unit width and a streamwise length l , 
the right hand side of (3.9) becomes, 

J 


h (T -T) = 
av aw ' 


I 


q d x 


(T - T) 
aw ' 


I 


d x 


r 


h d x 


(3.11) 
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as suming an isothermal wall (T f- T (x)). Incorporating (3. 10) in (3. 11) 
and carrying out the integration yields 

h = 0.664 Pr‘/> K 

av 

Using (3. 12) in (3.9) provides the average heat input which is equated 
to the net radiation heat loss, (2. 10), to obtain a fourth order algebraic 
equation for the isothermal equilibrium temperature on the flat plate. 

Laminar Flow Over a Cone 

Supersonic, conical flow is amenable to boundary layer calculations 
for the same reason as flat plate flow; that is, because along any ray 
the Mach number and the pressure are constant. Therefore, the 
preceeding analysis for the laminar flat plate can be used after account 
has been made for the boundary layer " stretching" as it proceeds down¬ 
stream to accommodate the ever increasing come diameter. The 
Mangier transformation can be applied to slender, axi-symmetric bodies 
to transform the three-dimensional flow into an equivalent flow in two 
dimensions. In Ref. 7 Covert shows that the Lighthill type of solution 
for conical flow is of the same form as that for the flat plate. Thus the 
equilibrium temperature distribution on a cone can be calculated using 
(3.5), (3.7) and (3. 8) ; however, in this case, 

(f-f 7F' 

It can be seen that the cone angle influences the temperature distribution 
only through setting the shock configuration, which in turn determines 
the local physical properties of the air. 

Now, to determine the average heat transfer to a laminar cone, 
the local heat transfer coefficient will still be based on the results of 
Ref. 6. The distortion of the boundary layer due to the conical nature 
of the flow can be accounted for by the Mangier transformation again. 

In Ref. 8 this transformation is applied, and the heat transfer parameters 
in the conical flow and equivalent flat plate flow cases can be related 


(-£H V 


(3. 12) 
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(3. 14) 


hflat plate 

The local heat transfer coefficient on the cone can be found by substituting 
(3. 10) into (3. 14) and evaluating the pertinent gas properties at the local 
free stream conditions. To determine an average heat transfer coef¬ 
ficient for the cone, an isothermal surface is assumed and the heat trans¬ 
ferred to the circumferential element of area at axial position, x, is 


written 


q dA = (T -T)'h2irsin6. xdx 
H aw 


(3. 15) 


where 6 is the cone half-angle. Integrating (3. 15) over the entire cone 
length, and dividing by the total surface area times (T aw “T) yields an 
average heat transfer coefficient; 


(0. 332) Pr 1 /* K 


(3. 16) 


This is the conical counterpart of (3. 12) and is used in the form of (3.9) 
to represent the average aerodynamic heat input to the cone. As in the 
exact solution, the cone angle effects the temperature only by adjusting 
the local free stream conditions. 


Turbulent Flow Over a Flat Plate 

The phenomenon of a turbulent boundary layer flow does not lend 
itself to neat solutions of relatively high accuracy because the mechanism 
of turbulence is not well understood, much less mathematically described 
adequately. Therefore, this discussion will be limited to a very simple 
approach to the turbulent boundary layer heat transfer problem with heavy 
reliance on experimentally observed trends and relations. Again the 
concept of a heat transfer coefficient and a temperature difference is 
used; and it is of use to introduce the non-dimensional Stanton number, 

St, defined as 
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(3.17) 


St = h/p u c 
' ~ 00 00 p 


where p w> and Cp are evaluated at local free stream conditions. 

(The subscripts used here are to avoid confusion with the same properties 
evaluated at the wall temperature which are of importance later.) 

In the absence of reliable heat transfer data over a large enough 
range of Mach number, the Reynolds analogy of a constant relationship 
between local Stanton number and local skin friction coefficient is used. 

In Ref. 9 it is pointed out that this relation is constant over a wide range 
of altitudes and Mach numbers and can be well approximated by: 


St = .6 c, 


Using this equation and the friction coefficient data presented in Ref. 9 
the ratio of compressible St to incompressible St is plotted in Fig. 2 as 
a function of Mach number. It is assumed that this curve is not changed 
significantly with variation in Reynolds number or heat transfer rate. 

Based on reasonably sound assumptions about the velocity distri¬ 
bution in the boundary layer, Ref. 10 gives the Stanton number for the 
incompressible, turbulent boundary layer as 

St. =(0.0288) (-E- jP~ ) 4/s Pr/* Kx‘‘ (p^Cpf* (3.18) 


To produce an isothermal approximation to the plate temperature as in 
the laminar case, the heat transfer rate is integrated over the length of 
the plate, then divided by the total area. Carrying out this operation on 
( 3. 18) yields , 


St. 


av 


0.036 Pr//;i Kl 

' II ' 


' 5 (p u c ) 
r oo oo p 


(3.19) 


Adjusting (3. 19) by the appropriate Stanton number ratio from Fig. 2, and 
evaluating p, p. and K at the wall temperature yields the average Stanton 
number and thus the average heat transfer coefficient to be used in (3.9) 
with the turbulent value of the adiabatic wall temperature. 
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Turbulent Flow Over a Cone 


1 


The Mangier transformation used in the discussion of the laminar 
cone can be modified to apply in the case of a turbulent boundary layer. 

In Ref. 8 it gives 


^cone 
S^flat plate 


1. 14 


(3. 20) 


as a result. By applying (3. 20) to (3. 18) the local Stanton number on the 
cone can be found. Using a procedure like that which produced (3. 16), 
the total heat transferred to the cone is found; assuming an isothermal 
condition, and then non-dimensionalized to obtain an average Stanton 
number. This process yields an incompressible equation, based on (3.18), 
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This expression is corrected for compressibility effects using the local 
Mach number and Fig. 2 to obtain the heat transfer coefficient for the 
isothermal approximation. It should be noted that, as in the laminar 
case, the cone angle controls the temperature only through adjusting 
the shock configuration and thus the local free stream conditions. 
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SECTION IV 

DISCUSSION OF APPROXIMATIONS AND EXAMPLES 

In this section the methods of the preceding two sections are 
applied to some particular cases to demonstrate the validity of the iso¬ 
thermal approximation, and to show the trends in radiation as some of 
the aerodynamic parameters are varied. The altitudes and velocities 
selected for the examples are somewhat arbitrary, however, it is felt 
that they cover the range of interest. A brief comment on the range of 
validity of the heat transfer equations seems to be in order here. The 
theories discussed in Section III are for perfect gas in continuum flow. 

At high Mach numbers the shock waves produced by blunt cones and plane 
surfaces at high angle of attack raise the air temperature to a level 
where the real gas effects must be considered. At gas temperatures 
above 2000°K the real gas effects become noticeable, and the methods 
of calculating heat transfer rates for those radiation calculations must 
be re-examined at temperatures above about 3500 °K. The assumptions 
involved in the continuum flow calculations begin to break down when 
M/4Re~ gets larger than about lO" 2 . It can be seen, however, that 
the range of validity of the methods discussed here covers Mach numbers 
up to seven and altitudes well above 100, 000 ft. ; certainly covering 
nearly all manned aircraft and nonballistic missiles. 

To evaluate the isothermal plate approximation in calculating the 
radiation from simple surfaces, a laminar, flat plate at zero angle of 
attack is considered. Assuming an emissivity of one and a length of 
three meters, the temperatures calculated by the Lighthill method are 
shown in Fig. 3 for Mach number s 3.0, 5.0 and 7.0, and at altitudes 
of 30,000 feet and 70,000 feet. The increasing equilibrium temper¬ 
atures with increasing Mach number is clearly evident. The 
lower aerodynamic heat transfer rates at higher altitudes (lower 
free stream density) show up as decreased equilibrium temperatures. 
The temperatures found by using the isothermal approximation are in¬ 
dicated on the curves. To get an accurate indication of the distribution 
and level of the infrared radiation from this plate, the Lighthill solutions 
for the temperature distributions were divided into a number of small 
isothermal sections, and the contributions of all of those, for a given 
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condition, were added together to obtain the total w^ vs. X curves. These 
"exact" results are shown, along with the results obtained by considering 
the plate to be isothermal, in Fig. 4 for the six conditions under consid¬ 
eration. 

These curves show the isothermal approximation to be very sat¬ 
isfactory at the low Mach numbers. At a Mach number of 3.0 the approx- 
imation leads to an error of less than 10% in w^ over the wavelengths of 
general interest. As the Mach number increases, the relatively high 
temperature region near the leading edge becomes more important and 
leads to some error in the approximation at short wavelengths. This 
can be seen in the curves for Mach 5. 0 on the short wavelength side of 
w \max‘ The total radiation predicted by the isothermal approxima¬ 
tion is about 90 percent of the more -exact value at the 30, 000 feet 
altitude, and only 80 percent of the exact value at 70,000 ft. At 
Mach number 7.0, even though the radiation level is well approximated 
at the long wavelengths (beyond about 4 microns), the short wavelength 
radiation is not well approximated by the isothermal solution, and the 
approximation predicts only about 80% andabout 66% of the total radiant 
energy at 30, 000 ft and 70, 000 ft respectively. This discrepancy is 
probably not great enough to warrant the additional labor involved in the 
use of the more exact temperature distribution, unless an accurate de¬ 
scription of the spectral distribution of the radiant energy is required 
for special purposes such as discrimination. 

To show directly the effect of changing Mach number and surface 
emissivity, the same flat plate is considered at an altitude of 30,000 ft. 
Figure 5 shows the results of calculating the radiation for three values of 
Mach number and three values of emissivity. As would be expected from 
the previous discussion, the energy level is raised, and \ max becomes 
smaller as the Mach number increases (due to an increase in heat transfer 
coefficient, and consequently a higher surface temperature). The families 
of curves showing the effects of changing emissivity reveal something 
rather unexpected — in general, the level of q\ max is not proportional 
t° « • At Mach 3.0 when the emissivity is reduced by a factor of two, 
q^ at the peak of the curve and in the long wavelengths is reduced by less 
than a factor of two while q^ below 3^ remains relatively unchanged. At 
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Mach numbers 5.0 and 7.0 the lp t o lOp radiation level is almost 
unchanged with changes in emissivity. This unusual effect exists because 
as the emissivity drops,the surface becomes a less efficient radiator, 
but the surface temperature rises. At higher values of e (lower surface 
temperatures) the aerodynamic heat transfer rate is higher [see (3. 9)] 
and consequently the rate of total radiation heat transfer is higher; how¬ 
ever, a greater fraction of this heat loss takes place in the longer wave¬ 
lengths and does not appreciably affect the region from lp to lOp. 

When the aerodynamic heat transfer rate is put into equation (2.9) 
to obtain an expression for the wall temperature, considerable simpli¬ 
fication can be obtained by omitting one or more of the radiation terms. 

In order to estimate the effects of these terms, the equilibrium tem¬ 
peratures ona 3 meter flat plate, having a laminar boundary layer and 
an emissivity of 0.7 were calculated for a range of Mach number and 
altitude. The results of these calculations are shown in Fig.6. In this 
figure the basic, equilibrium temperature solution, (T w ), against which 
the others are compared, includes only the radiation heat loss from the 
surface. The T' represents the equilibrium temperature found assuming 
other forms of the radiation portion of the heat balance equation. The 
total infrared radiation from the surface is proportional to the fourth 
power of the surface temperature, hence the (T'/ T w )^ vs. altitude 
presentation for various Mach numbers. 

If the cooling effects of radiation are neglected (3.9) equals 
zero and the equilibrium temperature equals the adiabatic wall temperatur 
Taw" ^ * s °t>vious from Fig. 6 that this cooling effect must be accounted 
for in all supersonic flight conditions. The ratio (T aw / T w )^ ranges 
from a modest 1. 3 to a gigantic 28, with the error becoming more 
pronounced as both the altitude and speed are increased. 

Two approximate curves showing the effects of heat radiated to 
the surface by the surrounding atmosphere also appear in the figure. The 
temperature increase over T w results in a total radiation increase of 
less than 10%atthe very high altitudes, and less than 5%.below 80, 000 ft. 
Clearly, this term in the heat balance equation need not be considered 
except in the most unusual situations. 
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Under some conditions the heat input due to the sun must be con¬ 
sidered. Using the radiation cooling and normal incidence solar radiation 
terms in the heat balance equations, the third family of curves in Fig. 6 
was produced. It can be seen that at low speeds and high altitudes the 
increased surface temperature due to the sunshine produces a significant 
increase in the total infrared radiation. 

As a final example, demonstrating the application of the methods 
discussed in this report, a typical aircraft shape will be considered. An 
airframe of approximately the size and proportions of a supersonic fighter 
is assumed and the radiation from the skin is calculated for two flight Mach 
numbers and two altitudes. Figure 7 shows the assumed configuration and 
it can be seen that the geometry is simply a combination of a cone, 
cylinders and rectangular flat plates for ease in calculation. It is felt that 
this simplification of geometry would not significantly affect the results 
of the analysis of most real aircraft. The surface emissivity is assumed 
to be 0.7 which is a representative value for most structural materials. 

In Fig. 8 the radiation is calculated for flight at Mach 2. 0 at 
40,000 ft. and three aspects, nose, plan and side. The contributions of 
the various structural components can be seen. In the nose aspect, the 
open area of the inlets is represented by a black surface at the free stream 
recovery temperature. The total target radiation curve is very near that 
of a gray body equal in area to the projected area of the target aircraft 
and at some near body temperature. This is true only at low Mach numbers 
where the isothermal approximation can be used. At high Mach numbers 
the high temperatures at short distance from leading edges produce large 
deviations in the radiation curves which are peculiar to each separate com¬ 
ponent and an over-all approximation can no longer be nade with any con¬ 
fidence. This phenonena is accentuated in the nose aspect, where the 
wing leading edges become very much hotter than any body mean tem¬ 
perature, and consequently add radiation at the shorter wavelengths. 

Figure 9 compares the radiation from the target aircraft at 40,000 ft. 
and 60,000 ft. and at Mach numbers of 2.0 and 3.0. Increasing the flight 
Mach number from 2. 0 to 3. 0 increases the total target radiation by about 
a factor of five. This is to be expected, of course, because of the in¬ 
creased adiabatic wall temperature. As the altitude is increased from 
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40, 000 ft to 60, 000 ft the total radiation level is decreased 20 % to 30.%. 
There is really no change in the adiabatic wall temperature, but the 
drop in free stream density causes a drop in heat transfer rate and there 
fore the radiation cooling becomes more effective. 
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Figure 5. Flux Per Unit Wavelength Versus Wavelength, 
Laminar Flat Plate at 30, 000 ft For Various 
Mach Numbers and Emissivities. f = 3 meters. 
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Lgure 9. Radiated Power Per Unit Wavelength and Solid 
Angle Versus Wavelength, Hypothetical Super- 


















APPENDIX A 

GRAPHICAL PRESENTATION OF PLANCK'S EQUATION 


It is desirable to have a simple, rapid method for graphically 
presenting the flux per unit wavelength as a function of wavelength for 
any grey body at any temperature. Equation (2.6) is too cumbersome to 
evaluate for each case and a method is outlined below which allows a plot 
of q\ versus X. to be made conveniently. 

Considering (2.6), unencumbered by the constants accounting 
for the geometry and emissivity, Planck's radiation equation can be 
written 

2hnc 2 ( ch/kXT , 

=—x?-V e - 1 


)- 


Now, if q^ and X are replaced by the variables m and n such that 


q K = Pm 
X = n/T 


(A-l) can be written 


m 


2hTTC* 

n* 



(A- 2) 


(A- 3) 


It can be seen that the variation of m with n is independent of temperature. 
A curve of m versus n, normalized so that the maximum value of m 
equals one, is shown in Fig. A-l. To obtain q^ and X the transformation 
(A-2) is used, i.e. the curve of Fig. A-l is multiplied by the appropriate 
values of T,m anda. (This just means changing the horizontal and ver¬ 
tical position of the curve on a log-log field without changing its shape) 

The relative values of any black body radiation curve can be 
obtained from Fig. A-l, and if one point on the curve can be determined 
absolutely the entire spectral distribution is known. The most convenient 
point to examine is the radiation "maximum" or the point of maximum 
q^. In Ref. 1 the wavelength at which this maximum occurs is given as 

X = he 

max (4.9651) kT (A-4) 
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or 


TX 

max 


2.898 X 10 3 ‘K(i 


and the value of at its maximum is given as 

q. = (3. 18 X 10 2 ) — ^J 5<r 
Xmax irhc 


(A-5) 


or 


Xmax = 1.2864 X 10* 


T* 


watts 
cm 2 °K^p. 


Thus, with equations (A-4) and (A-5) to determine the appropriate position 
on the q^“X plot of the radiation maximum, the curve of Fig. A-l can 
be laid over the log-log paper and the curve drawn for any black body. 

Of course the value of qXmax can be adjusted to account for emissivity 
and geometric considerations. 
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Figure A-l. Transformed Flux Per Unit Wavelength Versus 
Transformed Wavelength. 
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